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ABSTRACT 

In order to understand the processes that quench star formation within rich clus- 
ters, we construct a library of subhalo orbits drawn from ACDM cosmological N-body 
simulations of four rich clusters. The orbits are combined with models of star forma- 
tion followed by quenching in the cluster environment to predict colours and spectro- 
scopic line indices of satellite galaxies. Simple models with only halo mass-dependent 
quenching and without environmental (i.e. cluster-dependent) quenching fail to repro- 
duce the observed cluster-centric colour and absorption linestrcngth gradients. Models 
in which star formation is instantly quenched at the virial radius also fail to match the 
observations. Better matches to the data are achieved by more complicated bulge-disc 
models in which the bulge stellar populations depend only on the galaxy subhalo mass 
while the disc quenching depends on the cluster environment. In the most successful 
models quenching begins at pericentre, operating on an exponential timescale of 2 - 
3 Gyr, with the shorter timescale being a better match to disc colours as a function 
of cluster-centric radius and the longer being a better fit to the radial dependence 
of stellar absorption line indices. The models thus imply that the environments of 
rich clusters must impact star formation rates of infalling galaxies on relatively long 
timescales - several times longer than a typical halo spends within the virial radius of 
a cluster. This scenario favours gentler quenching mechanisms such as slow "strangu- 
lation" over more rapid ram-pressure stripping. 

Key words: galaxies: clusters: general - galaxies: formation - galaxies: evolution - 
galaxies: haloes - galaxies: star formation. 



1 INTRODUCTION 

Understanding the physical mechanisms that halt star for- 
mation has been a long-standing challenge in galaxy forma- 
tion and evolution. It is clear from the colour-magnitude 
relation l|Sandage fc Visvanathanl Il978l : iBower et all Il992l ) 

that mass or, more accurate ly, velocity dispersion 

i|Smith et al.l l2009al : iGraves et ail |2009h is a driving pa- 
rameter of the stellar populations in red galaxies. At the 
same time, the mo r phology-density dHubble fc Humasonl 
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The latter suggest that environ- 
star formation, moving galaxies from the 
blue cloud to the red sequence. Here quenching is a term 



meant to be evocative of or analogous to cooling a hot ob- 
ject. Quenching can be abrupt, as in a small campfire doused 
by buckets of water, or operate on longer timescales, as in 
cooling molten hot iron in a pool of water. It has become 
common to model both an "internal" quenching mechanism 
(assumed to be tied to galaxy mass) as well as environmen- 
tal quenching when a g alaxy falls into a higher mass halo 
and becomes a satellite dKauffmann et all 19931: Cole et all 



20001: Ivan den Bosch et all 120081; IWeinmann et al.l 1200a 
Peng et alj|20ld . l201ll : IWetzel et alj|2012bh . 



Many mechanisms that might quench star formation 
in cluster environments h ave been proposed: ram-pressure 



stripp ing of the cold gas l|Gunn fc Gott 1972 : Abadi et al.1 
Il999h . the removal of the hot gas halo ((Larson et al.lll98' 



some times known as "strangu lation" ([B alogh & Morri 
l2000h either via ram-pressure l|McCarthv et al.l 12008^ 
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via tidal strip ping by the cluster potential (|Merrittl Fl984; 
Mamonl 1 19871). or by "harassment" - e ncounters with othe r 
galaxies (fciallagher fe Ostrikeri Il972l ; iMoore et all Il998h . 
The ram-pressure stripping and strangulation mechanisms 
remove the gas that is the fuel for star formation, and so 
the stellar disc will dim and redden as its stars age. Sev- 
eral of the processes affecting l ate-ty pe galaxies have been 
reviewed bv lBoselli fc Gavazzil (|2006T ) . 

It has long been known that spir al galaxies in the cores 
of rich clusters are deficient in HI (IDavies fc Lewis! 11973: 
Havnes et al.lll984lGfovanerli fc Havneslll985l ; ISolanes et ail 



200 ll ) and stripping of star-forming gas has been observed in 
nearby rich clusters jKoopmann fc Kennev 2004] ; Sun et al 



2007] ; lYagi et al l 120071; 
ISive 



201C ; ISivanandam et ah! 



Yoshida et al.l 120081 ; ISmith et al . 
201Ch . usually appearing to strip 



gas on short (several hundred Myr) timescales, but it is un- 
clear if this is dominant quenching mechanism. Other, less 
abrupt or "gentler" mechanisms such as strangulation or 
tidal stripping are more difficult to observe. To determine 
which mechanisms are dominant, one might instead tackle 
the following more empirical questions: At what location in 
the cluster or group does quenching first occur? What is 
the timescale of quenching - is quenching abrupt, ending all 
star formation within tens of millions of years, or does star 
formation slowly decline over billions of years? What range 
of satellite galaxy masses are effectively quenched by the 
environment? Are all morphological components quenched 
equally, or, for example, are bulges insensitive to environ- 
ment while discs are quenched? 

A common approach to understanding th ese proc e sses is 
semi-analytic modelling (SAMs for short; see iBaughl (|2006T l 
for a review). SAMs attempt to model important physical 
process via simple analytical prescriptions, combining these 
with results from numerical experiments like N-body simula- 
tions (hence the semi). A disadvantage of the SAM approach 
when applied to all galaxy populations is that it is diffi- 
cult to isolate and study specific environments and physical 
processes, particularly with a large number of potentially 
correlated free parameters and physical processes. It is also 
challenging to determine if the parameterizations of physi- 
cal processes are themselves appropriate, since a wide range 
of values for the free parameters can allow for spurious, un- 
physical good fits to some observables. An alternative ap- 
proach is to focus on a particular environment with a more 
specific set of observables allowing for a limited number of 
mechanisms and prescriptions. 

In this paper, we will examine the link between star 
formation rates and the cluster environ ment. Previous mod- 
els of galaxy evolu tion in a clusters (|Balogh et al.l |2000| ; 
iDiaferio et alj|200ll and references therein) employed mod- 
els linking star formation rates with gas consumption , 
and in some cases, replenishment ([Weinmann et al,|[2 010). 
More recent studi es |Berrier et al.ll2009l ; lMcGee et al.ll2009l ; 
ISmith et al]|2012h have focused on pre-processing in galaxy 
groups. While these approaches have been successful in re- 
producing star formation rates (SFR) and colour gradients, 
they have not yet isolated the mechanism responsible for 
halting star formation despite including a wide variety of 
models for physical processes. We instead seek to produce 
models that employ an easily-tested causal connection be- 
tween environment (traced by infall into a r ich cluster) and 
the evolution of an individual galaxy (see iMahaian et al.l 



|201ll) for a recent example of star formation models of clus- 
ter galaxies). 

Rather than making an ab initio prediction of the total 
stellar mass content in each galaxy, as in most SAMs, we 
assign stellar masses by matching dark matter halo velocity 
dispersions to stellar velocity dispersions. The advantage of 
this approach is that the free parameters in our models can 
be given arbitrary values while still ensuring that each halo 
contains a galaxy of an appropriate stellar mass. We can thus 
assign arbitrary values to free parameters and indeed use 
any appropriate parameterization of star formation. How- 
ever, each model is not guaranteed to predict appropriate 
star formation histories for each galaxy, so our models are 
permitted to fail in reproducing some observations, unlike 
in most SAMs where the free parameters are tuned to re- 
produce many more observations. One drawback to this ap- 
proach is that it is difficult to model arbitrarily large num- 
bers of physical processes without an automated parameter 
search. This issue can be avoided by limiting the models to 
a very specific environment like rich galaxy clusters. If the 
stellar populations of cluster galaxies are deeply linked to 
their environment, one would expect star formation histories 
quenched by the cluster environment to successfully repro- 
duce cluster-centric radial colour and line index gradients. 
Conversely, if the cluster environment plays a limited role, 
models linking evolution to internal properties of a galaxy 
should more closely match observed data. 

In our models, galaxies are assumed to be forming stars 
until quenched by one of two mechanisms. In one scenario, 
"internal" quenching of star formation occurs at a time 
which is assumed to be a function of galaxy velocity dis- 
persion, such that galaxies with high velocity dispersions 
are internally truncated earlier (see Section H.l. II for details 
of our particular implementation). In the alternative sce- 
nario, environmental quenching of star formation occurs in 
the galaxy when its halo crosses some multiple of the virial 
radius of the cluster, with an option to delay quenching un- 
til the galaxy reaches pericentre or for some fixed time after 
infall. Though these implementations depict an overly sim- 
plified view of galaxy evolution, testing models based upon 
these mechanisms should provide clues as to the relative im- 
portance of environments in star formation. By varying the 
quenching radius of the cluster we also aim to test the range 
at which environmental quenching would need to be effective 
to reproduce observations. 

Most models to date have compared predictions with 
observations of galaxy luminosities and colours. However, 
the well-known degeneracies between age, metallicity and 
a-element enhancement can only be broken by using spec- 
troscopic absorption lines sensitive to these three parame- 
ters. Our aim is to break these degeneracies by comparing 
our data to observations of stellar absorption linestrengths 
in cluster galaxies. In this paper, we focus on colour and 
absorption lin estrength data from the NOAO Fundamental 
Plane Survey l|Smith et al.ll2004 hereafter NFPS) of galax- 
ies in rich clusters. 

This paper is structured as follows. Section[2]details the 
simulations used and presents information on the simulated 
clusters. Section [3] descibes the NFPS data in more detail. 
Section [4] describes the star formation models employed and 
the results of applying these models to galaxies. We discuss 
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the implications of our results in Section [5] and summarize 
and conclude in Section [6] 



2 N-BODY SIMULATIONS AND HALO 
ORBITS 

2.1 Methods 

The simulation data consists of four rich clusters ob- 
tained using the publicly a vailable parallel Tree-PM code 
GADGET-2 jSpringdl2()05h . Cluster candidates were iden- 
tified from a low-resolution initial simulation using 256 3 par- 
ticles in a cube of side length 512/i _1 Mpc. ACDM cosmo- 
logical parameters = 0.72, fi m — 0.28, h = 0.72 and 
as = 0.8 were adopted, con sistent with WMAP 7-year re- 
sults (|Komatsu et al.l 1201 ll ). Five candidate clusters were 
chosen for a higher-resolution re-simulation of 150 approx- 
imately equally-spaced time steps from z — 3 to z = 0. Of 
these five clusters, four were chosen for the final sample, 
with one candidate excluded due to a complex merger event 
between two massive clusters. The re maining four were re- 
simulated with the 'zoom' technique dKatz fc White! 1 19931 ; 
iNavarro fc White! Il994l ; IPower et all I2003T ). In brief, only 
particles passing near the central overdensity in each simu- 
lation are fully sampled. Particles which remain sufficiently 
far (in practice, more than about 20 Mpc/h) are subsampled 
and given higher masses to provide an accurate tidal field, 
with four levels of refinement up to the full 512ft _1 Mpc 
side box. Each re-simulation has 12.5 million dark parti- 
cles, 8.8 million of which are full resolution with a mass of 
6.16 x 10 8 M Q //i. The particle mass and minimum gravita- 
tional softening length of l/i -1 kpc are sufficient to resolve 
a Milky Way or M31 mass halo with at least 1,000 particles 
and achieve a resolution limit of 30 particles per halo for ha- 
los close to the total mass of a dwarf galaxy like the Large 
Magellanic Cloud. 

The masses and other physical properties of each cluster 
at z = are given in Table [T] Subhalo catalogues at each 
time step were generated using the AMIGA Halo Finder 
l|Knollmann fc Knebej2009l . hereafter AHF). AHF calculates 
isodensity contours on an adaptive mesh grid and identifies 
subhalos as collections of mutually bound particles within 
unambiguous density contours (i.e. overdensities) . The re- 
sulting subhalo statistics are for the simulated clusters are 
given in Table [1] AHF determines halo sizes from upturns in 
radial density contours, using an iterative process to remove 
unbound particles. We create halo merger trees by linking 
between consecutive snapshots' subhalo catalogues. 

A typical problem with N-body simulations of dissipa- 
tionless dark matter is the overmerging of low-mass dark 
matter subhalos as they are tidally disrupted by the clus- 
ter. This effect has b een shown to be at the least strongly 
resolution-dependent (|Moore et al.lll996 | ). and mo stly likely 
an artifact of poor resolution ( Klypin et al.|[l999l ). A com- 
mon solution to this "overmerging" problem, adopted here, 
is to use the most bound particle to track the orbits of 
these 'orphan' haloes l|Kauffmann et al1ll999l ; rSpringel et al.l 
l200ll ). In our implementation, any halo which is detected in 
at least one snapshot is tracked through its most bound par- 
ticles in all subsequent snapshots, even if it is not detected 
again by AHF. Thus once a halo is detected at high redshift, 
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a Millions of particles within the virial radius of the cluster 

^ Number of halos ever tracked, including recovered 'orphan' halos. 

c Number of halos which have ever crossed the virial radius of the 

cluster. 

^ Number of halos which crossed the virial radius but are now outside 
of it. 

Table 1. Properties of the four simulated rich clusters and their 
subhalo populations at z = 0. 



it is guaranteed to still be present in the final z=0 snapshot. 
Although AHF is capable of detecting most of the massive 
subhalos in the simulation, the most bound particle tracking 
technique recovers on average four times as many halos as 
in the original halo AHF catalogue. 

Although this method recovers many genuinely inde- 
pendent subhalos, a side effect is that subhalos are effec- 
tively prevented from merging. This is not an unreasonable 
assumption for most field galaxies falling directly in to a rich 
cluster. The large velocity dispersion in clusters will suppress 
mergers and the dynamical friction timescale is long enough 
that satellites can survive for billions o f years llTormen et al.l 
1998), unlike in lower mass clusters l|TormenT 19971 ). It is 
probably not a safe assumption for halos in smaller groups 
and likely results in an overestimate in the number of sub- 
halos near the cluster center, which could well have actually 
merged. Besides tracking overmerged halos, we use the most 
bound particles to fill in "gaps" in each halo's orbital his- 
tory, e.g. for halos detected in the 8th and 10th snapshot 
but not the 9th. Such gaps can occur when a halo crosses a 
particularly dense region of the cluster, temporarily evading 
detection before re-emerging in a less dense neighbourhood 
after pericentre. 

In this paper, we will typically define "infall" to occur 
when the subhalo crosses the radius r2oo, defined to be the 
radius at which the density exceeds the critical density, or 
in some models at a smaller fraction of T2oo- The lookback 
time at the first r2oo-crossing is noted as t r200 . The orbits 
allow us to also define the lookback time to first pericentre 
passage, t pcr . For the rich clusters studied here, we find that 



typically t pcl occurs approximately 0.5 Gyr after ^200 • tpcr 
grows roughly linearly with time as the cluster's dynamical 
time grows. Recently "infallen" halos take closer to 0.8 Gyr 
to reach pericentre on average, as the virial radius of typical 
clusters grows faster than the velocity dispersion. In princi- 
ple, pericentre can occur anywhere within r2oo; in practice, 
few halos pass closer than 0.1r2oo- The median pericentre 
is 0.4>2oo with no dependence on t in f a u except at very low 
redshift (Fig. [TJ, where it appears to dip to 0.3r2oo- 

In addition to tracking subhalo positions and velocities 
after infall, bound particles are used to compute velocity dis- 
persions for each subhalo. Subhalos are assigned a maximum 
velocity dispersion cr max corresponding to the velocity dis- 
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persion at the time when the subhalo had the largest mass 
(and in turn had the largest amount of bound particles) 
prior to infall within T2oo- This is based on the assumption 
that a subhalo will not grow significantly after infall, which 
is expected to be true. In the real universe, we might ex- 
pect the galaxy's baryonic component, which has dissipated 
and become more tightly-bound, to be much less sensitive 
to tidal disruption than a dissipationless dark matter halo. 
Thus the maximum pre-infall halo velocity dispersion serves 
as a proxy for the (post-infall) stellar velocity dispersion 
l|Conrov et alj|2006h . which in turn is the primary driver of 
the stellar populations (see Section U below). 

The simulations allow us to mimic observational sam- 
ples by viewing the subhalos positions and velocities in pro- 
jection. In all plots of projected quantities, we use the three 
axial projections of each cluster and appropriate line-of-sight 
velocities and velocity dispersions to average over projec- 
tion effects and improve galaxy statistics. We denote the 
projected radius on the sky as R, to distinguish it from the 
3D radius r. We also apply cuts in velocity-space, including 
all subhalos with velocity cz within ±3<t c i of the cluster's 
mean ~cz. Thus in the plots shown in projection, about 17 
percent of the plotted subhalos have R < 7-200 although in 
3D they are outside the virial radius (r > r2oo). This frac- 
tion can vary from 12-20 percent with velocity dispersion or 
infall-time related cuts, which we will apply when populating 
halos with galaxies in Section 2] For a more d etailed discus- 
sion o f contamination due to projection, see iMamon et al.l 
(|20ld ). 

2.2 Results 

We now examine cluster-centric trends of the two key 
subhalo properties to be passed to our galaxy evolution 
code in Section [4] subhalo velocity dispersion cr max and in- 
fall/pericentre time. In our models, these two quantities will 
determine the star formation history of the galaxy stellar 
populations. These quantities are entirely determined by the 
dynamics of the dark matter, as long as the assumption of 
a strict correlation between subhalo velocity dispersion and 
stellar velocity dispersion holds. Although this is challeng- 
ing to prove observationally, it is unlikely to be far from true 
given that most early-type galaxies are dynamically relaxed 
and hav e been found to ha ve nearly isothermal total density 
profiles (|Auger et al.ll2010fl . 

Fig.[T]shows the distribution of subhalo infall times past 
r2oo as a function of both projected and real distance to the 
centre of each cluster. Most recent infallen halos (< 2 Gyr) 
are clustered around the same position for any given infall 
time, likely because most halos have mostly radial orbits and 
thus have similar infall histories. The majority of intermedi- 
ate (2-4 Gyr) infallen halos are at large radii, even outside 
the virial radius of the cluster and thus form a "backsplash" 
population of galaxies that fell in and back out of the clus- 
ter. Few backsplash galaxies remain so longer than 6 Gyr; 
such old infallen halos are concentrated near the centre of 
the cluster. This virialized population of subhalos outnum- 
bers the very recent infalls within 0.25r2oo- At 0.5r2oo, the 
median infall time of a subhalo is barely 2 Gyr. 

When considering only infallen halos, the backsplash 
population dominates at large radii (> 0.5r2oo and the me- 
dian infall time increases once again. However, when includ- 



ing all halos (not just those that have fallen in), the median 
infall time decreases again past r2oo as the halo population 
is dominated by field galaxies falling in to the cluster for 
the first time, i.e. infalling halos. These halos can be as- 
signed zero or negative infall times based on estimates for 
their inevitable infall time onto the cluster (Fig. ??), where a 
negative median infall time in a particular bin indicates that 
more than half of the halos have yet to fall in. The tight lin- 
ear relation between these estimated infall times and r/r2oo 
suggests that most field galaxies within r2oo are falling in 
to the cluster on nearly direct radial orbits. As a result, 
despite the scatter in pericentres it is clear in 3D that the 
infallen halo population can be separated by radial velocity 
up to an infall time of 2 - 3 Gyr, where these recently in- 
fallen halos begin to reach apocenter. Halos with U200 > 4 
Gyr are almost all part of the virialized population, since 
few are backsplash galaxies and there is a roughly even split 
between approaching and receding galaxies. 

The trends identified in 3D are less evident in pro- 
jection. However, the clear distinction between approach- 
ing and receding halos remains. Similarly, there are two 
triangularly-shaped 'holes' in the projected phase space, 
with virtually no halos having infall times between 1-2 Gyr 
within 0.5i?2oo and very few backsplash halos with i r 2oo < 2 
Gyr. If there is a strong environmental impact on infalling 
galaxies, there should be two distinct populations past -R200 
- currently infalling halos hosting field galaxies and a smaller 
backsplash population hosting quenched galaxies. Within 
^200, the distinction is less clear - halos have a variety of 
infall times, and some are projection effects which have yet 
to fall in. These projection effects serve to lower the median 
infall time. As we will demonstrate in Section [4. II the trends 
identified in infall time can be mapped onto a stellar popu- 
lation age for halos hosting galaxies undergoing infall-based 
quenching, which is particularly appropriate for discs. In 
particular, the sharp drop in median infall time between 0.8- 
1 R200 will cause a similarly large change in optical colours 
of galaxies for models employing abrupt quenching on infall. 

Although there is considerable scatter, subhalo infall 
times are also correlated with the velocity dispersion and 
mass of the halo (Fig. The median infall time for dwarf 
halos (a < 50 km s _1 ) is 2-3 Gyr, while almost all massive 
subhalos (a > 150 km s" 1 ) fell in at least 2 Gyr ago with a 
median of 6 Gyr. This can be interpreted as a consequence 
of the spatial correlation/clustering of overdensities in the 
early universe, such that the largest overdensity collapsed 
into the cluster itself while nearby overdensities formed mas- 
sive subhalos which fell in at early times. The consequences 
of this correlation on galaxy formation models are that even 
a purely environmentally-quenched model would predict a 
non-zero slope in the colour-magnitude relation, since more 
massive (and hence brighter) galaxies would have fallen in 
earlier and hence be redder on average. 

In addition to the mass trend, there is a weak but non- 
zero trend in the median velocity dispersions of halos as a 
function of cluster-centric distance (Fig. |3J), even in projec- 
tion. This can be viewed as a combination of the infall time 
versus position (Fig. [TJ and infall time versus mass relation 
(Fig. 0, which essentially causes mass segregation of sub- 
halos. One would then expect a model without any explicit 
environmental dependence to maintain some age trend, as 
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Figure 1. Lookback time at infall t r 200 for subhalos as a function of radius, in three dimensions (left panel) and in projection (right 
panel). Halos are colour coded on whether they are approaching (blue) or receding from (red) the cluster. Halos that have not fallen in 
yet are given a naive estimate of their eventual infall time assuming no changes in their current radial velocity and the size of the cluster. 
Empty squares show the median infall times for all halos; squares with an x inside show median infall times for only those halos that 
have fallen in (t r 200 > 0). 




Figure 2. The median lookback times at infall of subhalos as a function of prc-infall velocity dispersion (left panel) and pre-infall mass 
(right panel). The squares show the median infall times. 
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R/R 



200 



Figure 3. The projected spatial distribution of subhalos as a 
function of velocity dispersion. The squares show the median ve- 
locity dispersion at various radii. Massive subhalos tend to be 
found near the cluster centre, where the median velocity disper- 
sion is ~ 20 km s — 1 greater than in the field. 



long as the quenching is tied to halo mass or velocity dis- 
persion. 



3 DATA 

Data for giant galaxies in clus ters are based on the NOAO 
Fundamental Plane Survey (ISmith et al ] |2004 hereafter 
NFPS) of ~ 3000 galaxies in 94 nearby rich clusters. Colours 
for 750 gia nt {Mr < —20 . 25) g alaxies in 8 rich clusters are 
taken from lHudson et all (|2010l ). This sample has no colour 
selection and includes red and blue galaxies. It provides (to- 
tal) galaxy colours in B-R as well as colours in the same 
filters for the bulge and disc components separately. Spectro- 
scopi c data, in the form of central absorption line indices, are 
from lNelan et all (|2005l ). The NFPS spectra are too noisy to 
obtain ages and metallicities for individual galaxies, but the 
trends as a function of internal velocity dispersion a and pro- 
j ected distance fro m the cluster centres are well determined 
l|Smith et al.ll2006l ). 



model stellar populations matching the observational lim- 
its of the data, i.e. Mr < —20.25. This typically results 
in a median velocity dispersion limit of > 100 km/s, or 
Miliar > MQ x 10 10 . The subhalos that host these galaxies 
have masses M ha io > 10 n M©. 



4.1 Single-component toy models 

We first explore two simple toy models designed to test the 
influence of cluster-centric environment. These are "single 
component" models in which the entire stellar population 
of the galaxy is described by a single star formation history 
(SFH). In Section 14.21 below we will explore more realistic 
bulge/disc models. 



4.1.1 The AHS Models 

We will use the lAllanson etafl (|2009l . hereafter AHS) SFH 
models as baselines for our single-component models. The 
AHS models parametrized the SFH in terms of simple func- 
tional forms (e.g. single burst, exponential, constant SF fol- 
lowed by "abrupt quenching", etc.) and then adjusted the 
parameters of the model in order to reproduce both the me- 
dian and the scatter in the global colours and central spec- 
troscopic line indices of NFPS r e d-sequence galaxies . Sev- 
eral studies l|Smith et al.l l2009bl ; iGraves et all 120091 1 have 
shown that stellar populations are more tightly linked to 
the velocity dispersion than to the stellar mass, so a priori, 
one expects a correlation between mean stellar age and the 
stellar velocity dispersion. In the AHS models, the central 
stellar velocity dispersion determines (1) an age-related pa- 
rameter, such as a Single Stellar Population (hereafter SSP) 
age or a quenching time, (2) the metallicity and (3) the 
a-enhancement. Random scatter is added to each of these 
parameters to generate a galaxy population with a distri- 
bution of ages, metallicities and a-enhancements. For each 
simulated galaxy, the stellar population parameters are used 
to generate simulated colours and absorption linestrength 
indices calculated b y convolving simple st ellar population 
(SSP) model SEDs llMarastonlll998l. boOFj) a nd Lick-index 
absorption linestrengths 1 Thomas et al.l 20031) with a given 
SFH in order to repr oduce the lin e indic es of red-sequence 
cluster galaxies from iNelan et all (|2005l ) and ISmith et al.l 
(2007). Having done so, AHS find reasonable agreement with 
galaxy colours as well. Furthermore, the AHS models include 
synthetic total magnitudes, allowing us to select samples in 
the same way as the real NFPS data. Thus by construction 
the AHS models match the observed red-sequence stellar 
population observables. 



4 STAR FORMATION MODELS 

We will now use the halo orbits of Section 12.21 as input 
to simple star formation models in which environmental, 
or more specifically cluster-centric "triggers" (such as in- 
fall past r2oo, or pericentre passage) quench star formation. 
This quenching can either be slow (e.g. causing an exponen- 
tial decay in the star formation rate), or abrupt. To gener- 
ate galaxy catalogues from the dark-matter only simulation, 
each subhalo is assumed to host a single galaxy, and the final 
simulated galaxy catalog contains only those galaxies with 



4-1.2 Internal Quenching versus Quenching on Cluster 
Infall 

The first ( "null hypothesis" ) model is one in which the clus- 
ter environment has no direct effect: we assume that the star 
formation history depends only on internal properties of the 
subhalo (specifically the velocity dispersion of the subhalo), 
and not on infall into the cluster. The subhalo's velocity 
dispersion before cluster-infall is used a proxy for the stellar 
velocity dispersion. Star formation rates are modelled either 
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as declining until a specified time tAQ or as a single instan- 
taneous burst (SSP). The a-enhancement and metallicities 
are assigned through a lookup table of similar galaxies in the 
AHS mock catalogues. This yields a red-sequence in which 
stellar population age, and hence galaxy colour, depends 
only the internal properties and not explicitly on position 
in the cluster. Nevertheless, colours of galaxies in this model 
have a weak and indirect dependence on cluster-centric ra- 
dius due to the correlation between subhalo velocity disper- 
sion and cluster-centric radius noted in Section [2] 

In the second simple model, the galaxy's internal ve- 
locity dispersion has no effect on its age and metallicity. 
Instead, the quenching time is set to the lookback time at 
which it fell into the halo, i.e. the time at which it first 
crossed T2oo- Metallicities and a-enhancements are set to 
solar with small (< 0.05 dex) scatter. To simulate this in 
the mock catalogues, we draw a galaxy from the AHS mock 
catalogues that has a similar tAQ as its infall time (indepen- 
dent of the galaxy velocity dispersion). 

Quenching can be accomplished in a variety of different 
ways. In these simple models, quenching is instantaneous 
and complete. Quenching can alternatively be modelled as 
an exponential decline in star formation rate with some post- 
infall e-folding time T post . Similarly, one might imagine that 
environmental effects are weak at T2oo and grow stronger 
at higher intra-cluster medium densities or simply higher 
densities for tidal stripping. To approximate this effect, we 
also consider models where quenching begins at pericentre, 
where densities should be highest. To briefly distinguish be- 
tween these models, Table [2] summarizes acronyms and gen- 
eral properties of these parameterizations, all of which will 
be elaborated upon. For the moment, however, it is instruc- 
tive to consider only the simplest mass-dependent MSSP 
and infall-dependent IAQ models to determine if further ad- 
ditions to the models are necessary. 

Fig.Ushows the comparison of the cluster-centric colour 
gradients for these simple single component models com- 
pared to the NFPS global, disc and bulge colours. The 
internally-quenched model has a very weak colour gradient, 
in contrast with the steeper dependence on cluster-centric 
radius in the total colours in the NFPS data. This weak 
gradient is, however, a much better match to the bulge 
colours, with the caveat that our models slightly underes- 
timate metallicities and produce an offset of about 0.05. 
The small variation in colours that is seen is the expected 
second-order effect produced by the correlation of position 
with velocity dispersion (Fig. [3|, in turn caused by mass 
segregation (see Section |2.2|) . If more massive galaxies are 
also redder, as they are in the NFPS data, a (weak) negative 
gradient with cluster-centric position results. 

A purely infall-based model will typically produce a 
wider range of galaxy colours and a steeper dependence of 
colour on cluster-centric radius. Galaxies that have recently 
(or never) fallen in have the bluest colours, dominating the 
outskirts of the cluster. In the core of the cluster, the time 
since infall spans a range from several hundred Myr for the 
most recently infallen halos still on their first orbit around 
the cluster core, to many Gyr for older infallen halos now 
bound to the cluster core, if not virialized. Hence the over- 
all colour gradient is a function of the relationship between 
infall time and position (Fig. [T]) modulated by the depen- 
dence of colour on infall time. The former relation is almost 



I 

m 




- MSSP 
IAQ 

- -x- - - PAQ 



0.2 0.4 0.6 0.8 1 

R A 2 00 

Figure 4. Single component model galaxy colours as a function 
of projected cluster-centric radius. These simple models include 
an old single stellar population, abrupt quenching on infall and 
abrupt quenching at pericenter. Black squares are NFPS total 
colours; blue are disc and red bulge colours. Error bars for the 
NFPS data show the uncertainty of the median colour of galaxies 
in each bin. 



entirely determined by cosmology, while the latter is depen- 
dent entirely on our choice of star formation model. These 
abruptly quenched models have steeper colour gradients as a 
function of cluster-centric radius than both total colours and 
disc colours, regardless of where quenching occurs (Fig. 

We conclude that these toy models bracket the range 
of possible formation histories for galaxies but are too sim- 
plistic to reproduce observed galaxy colours. We now turn 
to more realistic bulge/disc models that are more complex 
variants of these simple toy models. 



4.2 Bulge plus disc models 

4-2.1 Environmental Dependence of Bulges and discs 

The B — R colours of bulge and dis c components o f gian t 
galaxies in clusters were studied by iHudson et al.l l|201Cf ) . 
They found that bulge colours were redder than discs (as ex- 
pected) but that bulge colours were independent of cluster- 
centric radius. In contrast, the colours of disc galaxies de- 
pend significantly on cluster-centric radius and drive the 
weaker trend of global colours as reproduced above (Fig. Q. 
The lack of radial dependence of the bulge colours is there- 
fore similar to the predictions of the "internally-quenched" 
toy model above, whereas the radial dependence of the disc 
colours is similar to the infall-quenched toy model. This sug- 
gests that we build a composite bulge/disc model in which 
only the discs are affected by environment. 

In the NFPS data, the bulge-to-disc ratio varies along 
the red sequence, with higher-mass systems being more 
bulge-dominated. Indeed AHS showed that the bulge-to- 
total light ratio (B/T) was tightly correlated with central 
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Model Name 


Acronym 


Quenching Method 


Quenching Model 


Quenching Position 


Pre-Infall r 
Gyr 


Post-Infall r 
Gyr 


Simple Stellar Population 


MSSP 


Mass 


SSP 


N/A 


N/A 


N/A 


Quenched at Infall 


IAQ 


Infall 


AQ 


r/r 20 o = 0.125 - 1 


oo — 1 





Quenched at Pericentre 


PAQ 


Pericentre 


AQ 


r/r 200 < 0.125 - 1 


CO — 1 





Strangled at Infall 


IS 


Infall 


STR 


r/r 200 = 0.125 - 1 


CO — 1 


0.5-7 


Strangled at Pericentre 


PS 


Pericentre 


STR 


r/r 2 oo < 0.125 - 1 


CO — 1 


0.5-7 



Table 2. Basic Star Formation Model Parameters. 



velocity dispersion (see their Figure Al). We use their empir- 
ical relation to determine B/T at a given velocity dispersion 

(7. 



4.2.2 Bulge Model 

Having fixed the bulge fraction, bulge ages are then gener- 
ated as for the "internally-quenched" toy model described 
above. The bulge metallicity is determined by choosing the 
total metallicity as a function of its velocity dispersion, by 
drawing galaxies at random from the AHS models. In prac- 
tice this means that the bulge has super-solar metallicity and 
is increasingly metal-rich at larger velocity dispersions. By 
contrast, discs are solar metallicity with small (< 0.05 dex) 
scatter and no a-enhancement. Since the total metallicity is 
a mass-weighted average of the disc and bulge metallicities, 
we adjust the bulge metallicity to reproduce the desired to- 
tal metallicity, with a corresponding increase in the bulge 
ct- enhancement . 



4-2.3 Disc Quenching Model 

The above prescriptions fix the stellar populations of the 
bulge and also the metallicity and a-enhancement of the 
disc. This leaves only the SFH of the disc to be mod- 
elled. We first consider disc colours in the field, where there 
is assumed to be no quenching. For field discs, we adopt 
an exponentially-decaying SFH and constrain the e-folding 
timescale so that the colours match t he median colour o f 
discs in the field: B-R = 1.25 ± 0.05 (|Hudson et al.ll2010D . 
This is in reasonable ag reement with the colours found by 
iMacArthur et all l|2004h for discs in field spiral galaxies in 
the same magnitude range as the NFPS sample. We find 
that the observed colours can be fit with a (pre-quenching) 
exponential timescale r prc = 4-5 Gyr which yields B — R 
in the range 1.28 - 1.20. 

In our bulge/disc models, disc quenching will be caused 
by interaction with the cluster environment. There are two 
parameters that control the quenching. The first parameter 
controls how rapid the quenching is. In the simplest scenario, 
star formation is instantly and completely suppressed, as 
might be expected if the quenching is due to ram-pressure 
stripping of cold gas. We refer to this as "abrupt" quenching 
(AQ). Alternatively, star formation may be quenched more 
gradually, as in the "strangulation" scenario. We model this 
as a second exponential decay in the star formation rate. The 
exponential decay time after quenching, r post , is uncertain 



and so below we experiment with different values. Note that 
in this context AQ corresponds to r post = 0. 

Two additional parameters control how quenching pro- 
ceeds. One parameter determines where in the cluster 
quenching first occurs: at r 2 oo (I) or at pericentre (P). An 
additional parameter "R" determines how far the influence 
of the cluster stretches as a factor of r 2 oo- By default and in 
most models this parameter has a value of one; if not oth- 
erwise specified this default value applies. In some models, 
however, R can be smaller than one. In the case of R — 0.5, 
only halos which cross 0.5r2oo are considered as infallen and 
hence quenched, so halos with pericentres at 0.6r2oo would 
not be quenched. 

4-2.4 Abrupt Quenching or Strangulation? 

We first examine the simplest case of abrupt quenching 
(AQ), setting the pre- infall e-folding time to r prc = 4 Gyr. 
For the moment, we compare only pericentre-quenched mod- 
els; we will examine the differences between pericentre and 
infall quenching shortly. We find that the abrupt quenching 
at r 2 oo produces a very sharp transition from blue (star- 
forming) to red (quenched) galaxies just within the quench- 
ing radius, as shown in Fig. [5] We have experimented with 
different quenching radii (e.g. at a fixed fraction of r^oo, or 
at pericentre) and find that changing the quenching radius 
does not alter sharpness of the transition from blue median 
colours to red, it merely shifts radius at which the sharp 
transition occurs. Such a sharp transition is not observed in 
the NFPS disc colour data. Thus, it is not possible to match 
the distribution of colours in the colour versus projected po- 
sition plot with any model that quenches star formation in- 
stantly, regardless of where the quenching takes place. This 
suggests that the disc star formation is quenched more grad- 
ually, as might be expected in a "strangulation" model. Fur- 
thermore, the location of the transition in median colours is 
not necessarily the same as that of the quenching radius. 
Since most galaxies reach pericentre at 0.4r2oo, backsplash 
galaxies can dominate out past 0.5r2oo- Only by shrinking 
the quenching radius below 0.4r2oo can a significant change 
in the transition point be achieved. 

Fig. [6] compares strangulation models in which the 
disc quenching occurs at pericentre (assuming pericentre 
is within r2oo) and with different values of pre- and post- 
quenching r. These models produce roughly linear median 
disc colour-position relations, bracketing the slope of the ob- 
servations. Note that the choice of r pro fixes the median disc 
colour of the outermost bin of field-dominated galaxies while 
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Figure 5. The left panel shows median disc colour as a function of projected cluster-centric radius for the NFPS observations (black 
squares as in F'ig. [4jl and "abrupt" (T pos t = 0) quenching models. The right hand panel shows the distribution of disc colours for quenching 
at 0.5r2oo illustrating the sharp transition from star forming to quenched galaxies. Blue points represent galaxies approaching the cluster 
center in 3D while red are receding. 



only slightly altering the colour of the innermost, quenched 
galaxy-dominated bin. In other words, central cluster galax- 
ies retain little memory of their initial conditions in these 
models; instead, their colours are most sensitive to r pos t, the 
speed of the quenching after infall. Models with r pos t < 2 
Gyr produce discs that are too red compared to galaxies in 
the innermost radial bin. The model with the closest match 
to the observations has T pos t = 2.5 Gyr, only slightly faster 
than the pre-quenching r prc = 4 Gyr. This suggests that on 
average quenching need not be especially strong to repro- 
duce the observed trends, and it certainly does not appear 
to be "abrupt". However, some quenching must be applied, 
as models with identical r prc and r post produce no colour 
gradient whatsoever. 

The strangulation models also a reasonable job repro- 
ducing the entire disc colour distribution and do not show 
the strong bimodality in galaxy colours exhibited by the 
abrupt quenching models in Fig. [5] Such a bimodality could 
be accommodated in a model with larger dispersions in 
galaxy metallicites and/or r prc , which might also be neces- 
sary to produce some of the very blue observed discs. Sim- 
ilarly, the observational unc ertainties on individ ual galaxy 
disc colours of ~ 0.1 mag |Hudson et al.l [20101 ) precludes 
the rejection of any model with strong colour bimodality. 
However, strangulation models are particularly attractive 
because the bimodality visible in colours as a function of 
3D cluster-centric distance is virtually completely erased in 
2D distances by projection effects, and so no additional ad- 
justments are required to produce a relatively smooth ra- 
dial dependence of galaxy colours on projected distance. 
Note, however, that independent evidence exists for strong 
bimodality in specific star formation rates of galaxies, which 
we will address in Section f4. 2. 61 



4-2.5 Quenching Location in the Cluster 



In the models discussed above and presented in Fig. [6] the 
strangulation is assumed to occur at pericentre, where the 
cluster's tidal field (and tidal stripping) is strongest. It is 
possible that a different mechanism is at play. For exam- 
ple, ram pressure stripping is sensitive to the density of 
the intracluster medium, which depends primarily on radius. 
We therefore consider whether strangulation models could 
be compatible with infall-based quenching at radii smaller 
than r2oo- Figure [7] shows the median colours for the best 
fit with T( pre , p0 at) = (4, 2.5) model. Strangulation at 0.5r2oo 
produces virtually indistinguishable results from strangula- 
tion at pericentre. Both models predict mostly blue discs 
near the outskirts and red discs near the cluster center. The 
model with strangulation at 0.5r2oo produces slightly bluer 
colours between 0.3 — 0.4r2oo, where halos are more likely 
to have fallen in to 0.5r2oo than to have reached pericentre, 
and so are slightly younger with an infall quenched model. 



Infall-quenched models with very small quenching radii 
(r q uench < 0.5r2oo) tend to produce a steep transition from 
the red cluster core to blue, unquenched galaxies. The me- 
dian colour becomes fiat at about 1.5 times the quenching 
radius. We note that the observations are not inconsistent 
with such a flat slope from 0.6 r2oo out to the field, in part 
due to poor statistics on the field population and large error 
bars on their median disc colours. Thus, we can only safely 
exclude models with quenching radii much smaller than 0.5 

^200- 
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Figure 6. The left panel is as in Fig.[5]but the models are now "strangulation" models in which the disc is exponentially quenched after 
reaching pericentre. The right-hand panel shows the distribution of disc colours for the PS-T(4,2.5) model. The post-quenching "tail" of 
star formation (strangulation) results in a steady decrease in median colours of galaxies. 
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Figure 7. Left panel as in Fig. [5] but now for "strangulation" models invoked at different quenching radii. Only galaxies falling within 
some fraction of the virial radius are quenched; the smaller the fraction, the closer to the cluster centre the transition from blue- to 
red-dominated galaxies occurs. The middle and right-hand panels show the distribution of galaxy colours for the IS-T(4,2.5)-R0.5 model 
as a function of radius (projected and 3D, respectively), illustrating how projection effects wash out the sharp transition in galaxy colours 
visible in three dimensions. 



4-2.6 Delayed Quenching and Bimodality 

In a similar study to this one, IWetzel et al] (|2012bh found 
good agreement with SDSS observations of specific star 
formation rates (SSFRs) using models with "delayed-then- 
rapid" quenching. The rapid quenc hing preserve s the b i- 
modality in the SSFRs observed by IWetzel et all l|2012al ). 
whereas the delay is necessary to match the overall quenched 
fraction as a function of environment. Their best fit model 
has a quenching that begins ~ 3 Gyr after infall, followed by 
rapid quenching with r post ~ 0.5 — 1 Gyr. We implement this 
by adding a constant delay time to the quenching time. The 



effect is similar to pericentre-based quenching, but with a 
longer delay after infall: pericentre occurs after at most ~ 1 
Gyr. After 3 Gyr, galaxies are in the "backsplash" popula- 
tion. 

Fig- H] compares the "delayed-then-rapid" quenching to 
the strangulation models. The two 3 Gyr delay models (la- 
belled with "-D3") wit h r pos t = 0.5 — 1 ro ughly bracket 
the range considered bv lWetzel etafl (|2012bl) . As expected 
the additional delay keeps more galaxies on the "blue se- 
quence" compared to a strangulation model with the same 
Tp OSt . A value of T post ~ 1 washed out the abrupt transition 
of the median colour from blue to red. But a weakness of 
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Figure 8. The left panel is as in Fig.[5]but now for "strangulation" models with quenching delayed for a time D after infall/pcricentre. 
The right-hand panel shows the distribution of disc colours for the IS-T(4,1)-D3 model, which displays qualitatively similar behaviour 
to Fig. [6] but with stronger bimodality and a redder red sequence. 



these models is that the mean slope of the B — R vs R/r2oo 
is ~ 0.2 magnitu des, whereas the dat a favour a shallower 
slope 0.1 ± 0.025 (|Hudson et alJl201Ch . This steep slope is 
found for any model with T pos t < 1, whether delayed or not. 

It is interesting to note that the general behaviour of 
delayed-then-rapid quenching models is to produce a steep 
gradient in median colours close to the cluster center and a 
flat slope at larger distances. This is similar to the behaviour 
found in quenching at small radii (Fig. [7}. One might then 
expect to be able to produce similar results in a model with a 
small quenching radius, smaller r pos t and a delay time in the 
form of pericentric quenching rather than infall quenching. 
In fact, the first such model tested produces very similar 
results to the delayed-then-rapid quenching models, as show 
in Fig. [8] While these small- radius quenching models do not 
produce much better fits to disk colours, they do show that 
one can achieve similar colour bimodality without invoking 
an arbitrarily long delay in quenching after infall. 

In addition to median colours in each bin, we com- 
pare the fraction of blue galaxies in different disc models 
in Fig. 1101 Although the choice of division between blue 
and red is somewhat arbitrary, the blue fractions are similar 
for all of the models previously shown in Fig. [8] increasing 
from 15-20 percent near the cluster core to 70-80 percent at 
7?2oo- The choice of colour cut (B-R of 1.32) can alter the 
fractions but produces similar results as long as it is between 
blue sequence and the red galaxies near -R200 , roughly divid- 
ing star-forming from quenched galaxies. The two models 
with the most similar median colours have virtually iden- 
tical blue fractions (PS-T(4,1.5)-R0.5 and IS-T(4,0.5)-D3), 
but curiously two models with quite different trends in me- 
dian colours also show very similar variations in blue frac- 
tion (our fiducial PS-T(4,2.5) model and IS-T(4,1)-D3). In 
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Figure 10. Blue fractions for the disc models in Fig. [9] defined 
as the fraction of galaxies with colours bluer than the cut of 1.32 
separating the blue from the red sequence. 



principle one could use the radial trends of blue fractions as 
an additional diagnostic to validate models. In practice, the 
uncertainties in observed disc colours are much too large for 
this exercise, which would also require more careful mod- 
elling of the intrinsic scatter in disc colours. 
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Figure 9. The left panel is as in Fig.[8]but now with additional "strangulation" models with a shorter T pt) st and a small cluster quenching 
radius. The right-hand panel shows the distribution of disc colours for the PS-T(4,1.5)-R0.5 model, which displays very similar behaviour 
to the right panel of Fig. [8] 
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Figure 11. Colours of bulge+disc models compared with NFPS 
total colours (black points). 



4-2.7 Consistency check: total colours 

We have shown that a viable model that reproduces both 
the bulge and disc colours, and their radial dependences, 
is one in which the bulges are quenched "internally" while 
discs are quenched by the cluster environment either at peri- 
centre or around 0.5 r2oo and on a timescale of T pos t ~ 2 - 3 



Gyr. We can make a further consistency check of this model 
by comparing total colours of simulated galaxies with NFPS 
rather than examining discs and bulges separately. Fig. QT] 
shows the median-colour - projected-distance relation for 
various two-component models. While none of the bulge- 
plus-disc models identically reproduce the observed data, 
they do roughly bracket the range of possibilities, given that 
the data are not dust-corrected and likely intrinsically some- 
what bluer than shown. 



4.3 The effect of dust 

Gala xy colours ar e , of c ourse, sensitive to the presence of 
dust. iDriver et al, I (|2007h find that, on average, bulge+disc 
galaxies in the field have a some attenuation of B-band light, 
well-descri bed by a central face-on opacity ~ 4 in the 
models of Tuffs et al.l l|2004l ) . It is we ll known that spira l 
galaxies in clusters are stripped of HI (|Solanes et al.ll200ll 1. 
If dust is also stripped, then this would introduce a differ- 
ence in colour solely due to the absence of dust in the clusters 
co mpared to i t s pres ence in field galaxies. Using the models 
of lTuffs et ail (2004), we estimate that if all the dust was 
completely stripped, the effect on disc colour would be to 
make the discs bluer by 0.10 to 0.22 magnitudes in B — R 
(with a median of 0.14 for a typical disc with an inclination 
i — 60°). Moreover, if the discs on the outskirts were not 



stripped (and hence contained dust and were redder) and the 
ones in the centre were stripped (and hence bluer), then the 
effect would be to make observed gradients shallower. Thus 
it is possible that the intrinsic (dust-free) cluster-centric gra- 
dient in disc colour is considerably steeper than the observed 
one. 

There are, however, two arguments against such an ex- 
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treme scenario. First, according to the models of lTuffs et al.l 
(2004), the bulge light (and bulge colour) are also strongly 
affected by dust. If the dust is stripped then the bulge will 
also appear bluer, and moreso for bulges than for discs. In 
their models, for a typical galaxy with a disc inclination of 
i = 60° the bulge is reddened by 0.35 mag. Thus, if strip- 
ping were important, we would expect bulges in the cluster 
cor e to be bluer than tho se at the cluster edges by 0.35 mag. 
Yet iHudson et al l (|2010h show that there is no gradient in 
bulge colour as a function of cluster-centric radius. While it 
remains possible in principle that an intrinsic dependence of 
bulge stellar populations exactly cancels the putative "dust 
removal" gradient, this seems rather contrived. Rather it 
suggests that the dust is not stripped in the cluster core. 

S econd, it is also pos sible to study dust stripping di- 
rectly. ICortese et alj (|2010T l studied dust stripping in Virgo. 
They found that there was evidence of dust stripping from 
the outer edges of discs, but the amount of dust stripping 
was high only in the most extreme "Hi-deficient" galaxies 
(defHI > 1). Galaxies with such high levels of Hi-deficiency 
are very rare in rich clusters, comp rising only a small frac- 
tion of the spiral galaxy population (jSolanes et alj200lh . We 
conclude that there is no strong evidence for dust stripping 
due to the cluster environment in the bulk of the infalling 
spiral population. 



4.4 Comparison with stellar absorption line 
indices 

Stellar absorption lines give detailed information on the ages 
and metallicities of non-star-forming stellar populations. In 
this section, we will compare the models to the NFPS data 
for red galaxies. Note that the NFPS spectra are fibre- 
spectra of the central regions (2 arcsecond diameter fibres). 
We can predict how much of the light within this central 
aperture arises from the bulge, and how much is from the 
disc, by usin g our bulge/di s c mod el. This calculation has 
been done by Hudson et alj (|2010l . see their Appendix Al) 
who derive an "aperture" B/T as a function of a, which 
ranges from 75-100 percent (depending on velocity disper- 
sion) for the NFPS sample. We adopt their empirical scaling 
here. Given the fraction of light from the bulge and from the 
disc, we simulate the central spectral line indices observed 
by NFPS. 

The NFPS line index data are for red galaxies only. 
To exclude the blue cloud in our models, we reject model 
galaxies bluer than 0.25 magnitudes from the red sequence, 
which is assumed to have a slope of 0.05/mag. In most cases, 
the presence of a significant fraction of bulge light ensures 
that very few galaxies are bluer than this cut. 

The median line-index-projected position relation is 
shown in Fig. 1121 H/3 line index trends are consistent with 
observations for the bulge-plus-disc strangulation models. 
The H,3 data are rather noisy, but the model that fits best 
is the one longest with the longest r post ~ 3 Gyr. The H7 
(HgF) data appears to disfavour many of the models. Only 
(pre, post) = (4, 3) model appears to predict 



the intermediate r. 



In summary, the line indices favour slightly longer post 
infall strangulation times than do the colours. Recall, how- 
ever, that the fibre data measure only the central 1.9 kpc. 
We have extrapolated the disc light assuming that it fol- 
lows an exponential profile right to the centre of the galaxy. 
The assumption may not be valid, particularly for SOs and 
early-type spirals. Furthermore, the Balmer line indices are 
sensitive to recent star formation, so what is relevant here 
is actually the scale lengths of the star-forming part of the 
exponential disc, which may differ from the fl-band scale 
length used in the bulge-disc model. iKoopmann et alj (|2006l ) 
found that Ha was more extended than i?-band light, at 
least in late-type galaxies in the field. Similarly, UV and 
blue colour gradients also suggest that, on average, recent 
s tar-formation is occurring preferentially in the outer disc 
Munoz-Mateos et al.ll2007h . particularly in Hi-rich galaxies 
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a sufficiently shallow radial gradient to at least marginally 
agree with the H7 data. In particular, the models with short 
quenching times (r pos t ~ 1 Gyr) predict very steep H7 gra- 
dients, as the virialized population will have virtually no 
young stars to produce strong H7 absorption. 



The situation for earlier-type galaxies 
and cluster galaxies is the less clear, but does suggest that 
the possibility that the central regions have less star-forming 
disc light than the simple exponential extrapolation used 
here. 



5 DISCUSSION 

We have shown that models with relatively rapid quenching 
fail to match the observations. By relatively rapid quench- 
ing, we mean either "abrupt" quenching, which corresponds 
to Tp OS t = 0, or "strangulation" models with r post < 1 Gyr. 
Such models fail in two ways. First, for galaxies in the cen- 
tres of clusters, these models produce discs which are redder 
than observed, and a cluster-centric radial gradient in disc 
colours which is too steep (Fig. [5]). Second, they overpredict 
the Balmer lines (Fig. I12[) . These basic trends are present 
whether the quenching occurs at infall or a fraction of r2oo 
or is delayed. 

Models in which star formation is quenched more 
slowly (2-3 Gyr) do a better job at reproducing both the 
disc colour cluster-centric radial gradients and the absorp- 
tion line trends. The colour data prefer timescales at the 
shorter end of the range (~ 2 Gyr) whereas the absorption 
linestrengths prefer > 3 Gyr timescales. 

It is more difficult to determine where in the cluster the 
quenching is triggered. Models in which strangulation begins 
at the virial radius or at pericentre are equally good fits to 
the disc colour gradients, but the latter has better physical 
motivation. However, the difference between otherwise iden- 
tical models employing quenching on infall vs pericentre is 
typically quite small, and the models can easily accomodate 
quenching at some intermediate point in the orbit as well. 

We can also compare the quenching times in our best- 
fit models to those quoted in other recent work s em- 
ploying environmental quenching. IWang et all (|2007tl mod- 
elled quenching of satellite galaxies, favouring an e-folding 
time (i.e. r DOB t of 2.5 Gyr, similar to our own best fits. 
|Pe Lucia etail (|2011r i presented models suggesting that 
satellites on average spend 5-7 Gyr within a halo of mass 
m > 10 13 Mq before being quenched. Although such models 
are not directly comparable to our own, as we only allow the 
cluster itself to quench subhalos, a 5-7 Gyr quenching time 
is not entirely inconsistent with some models depending on 
when one considers a galaxy to be quenched. For example, 
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Figure 12. Differential stellar absorption line indices as a function of projected position. The left panel shows H/3 and the right panel 
shows H -v. The vertical ax is is the residual from the mean line index-cr relation, normalized to at R = 0. NFPS data for giant red 
galaxies llSmith et al. 2006) arc shown by the black squares with error bars. The curves show bulge-plus-disc models for red galaxies with 
disc strangulation at pericentre, but differing pre- and post-infall timescales. For example, the black solid curve labelled M+-PS-T(4,3) 
has a r pre = 4 Gyr and T post = 3 Gyr and is a good fit to the H/3 cluster-centric radial gradient, and only slightly overpredicts the H7 
gradient. The other curves have shorter post-infall quenching times and are poorer fits. 



in the T post = 3 Gyr model, it would take nearly 7 Gyr for 
star formation to drop to 10 percent of its pre-quenching 
value after reaching pericentre. 

Several recent studies (|Weinmann et all l201fj : 
IWetzel et al.l l2012bl ) have attempted to reproduce the 
observed bimodality in specific star formation rates (SS- 
FRs) in SDSS galaxies with (respectively) full semi-analytic 
models or simple parametric star formation models. The 
generally-favoured model with r post ~ 2.5 Gyr does not 
produce a strong bimodality in disc or global galaxy colours 
as a result of the rel atively long e - folding time. For this 
reason, the models of IWetzel et al.l l|2012bh prefer a short 
quenching e-folding time, equivalent to our T pos t. To avoid 
producing too many red galaxies, they introduce a ~ 3 
Gyr delay before quenching begins. However, as noted 
above, Fig. [8] shows that the short r post also produces discs 
in cluster centres that are too red compared to those at 
the virial radius. They also predict H7 linestrengths for 
red galaxies near the virial radius that are stronger than 
observed (Fig. [12}. 

Although most of our data disfavour models with short 
r pos t and long delay times, we have been able to pro- 
duce alternative models with similarly (but slightly weaker) 
colour bimodality (Fig. [9}. This is achieved by shrinking 
the quenching radius, quenching at infall instead of pericen- 
tre (producing a short delay), and increasing T pos t slightly. 
These models are perhaps more physically plausible as they 
do not introduce an arbitrary delay time and there are no 
particularly strong constraints on whether the quenching 
radius is closer to 0.5 or 1 x r2oo- Still, these alternative 
models do not provide much better fits to the colour or line- 



index trends, so they do not solve the general problem that 
the relatively rapid quenching necessary to produce bimodal 
colours produces discs that are too red. 

There are several ways in which these studies could 
be improved. The overmerging problem can be addressed 
partially with better numerical resolution. More recent halo 
finding algorithms use partic le velocities in add ition to po- 
sitions and binding energies (|Oman et al.ll2012r i. which can 
improve estimates of infall and pericentre times and pericen- 
tre distances. Both of these improvements would likely affect 
low mass halos containing dwarf galaxies more so than the 
massive halos in our galaxy sample, which are well resolved. 
Halo velocities can also be use d to disentangle proje ction 
effects and backsplash galaxies (|Mahaian et al.l (|2011h ). al- 
though large observational samples are required to detect 
significant differences in the data. 

Because of our focus on reproducing mean/median 
trends, we did not include many ingredients into our disc 
models, such as mass-dependent SFRs, and scatter in 
the SFHs which one might model as variable r pre . These 
ingredients would likely help in reproducing the scatter 
about the mean relations and better reproducing mass- 
dependent trends, e.g. for dwarfs vs. giants, for which 
we do note substantial differences in the infall history. 
Another missing ingredient in the models is variable or 
density-dependent quenching strengths. Quenching mech- 
anisms like ram-pressure stripping can depend strongly 
on intr a-cluster gas density, pressure and satell ite infall 
velocitv (|Gunn fc Gott]|l972l ; iMcCarthv et al. 2008). The in- 
clination of the satellite's disc may also have an impact. 
These effects could introduce significant scatter in the val- 
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ues of Tpost and explain the more extreme observations of 
ram-pressure stripping in galaxies. 

The models presented here only consider quenching 
from the cluster itself and neglect pre-processing in smaller 
groups prior to infall. However, Fig. [2] shows that subhalos 
massive enough (a > 200km /s) to pre-process their own 
satellites fell in over 5 Gyr ago. Thus, the galaxies most 
likely to have been affected by pre-processing are already 
part of the old, virialized population. Furthermore, our best- 
fit disc models are tuned to reproduce field galaxies outside 
of galaxy clusters. These field galaxies include (potentially) 
pre-processed disc galaxies in groups, and so the disc com- 
ponents in our models reproduce the mean effect of pre- 
processing by construction. 

Despite these caveats, it is rather remarkable that rea- 
sonable agreement is found between different models and 
observational data sets. It is also surprising that models 
with as few as two free parameters can match galaxy colours 
and absorption linestrengths from cluster centre to outskirts, 
though simultaneously matching both observations is diffi- 
cult. All of our successful models favour relatively gentle en- 
vironmental quenching with timescales longer than the dy- 
namical time in the cluster and longer than typically quoted 
for ram-pressure stripping of cold gas discs. 



6 SUMMARY 

We have constructed a catalog of subhalos in dark matter 
simulations of four rich cluster of galaxies. By tracking halo 
orbits, we have demonstrated that the infall history of halos 
imprints a relationship between two key parameters (infall 
time and halo velocity dispersion) of halos and their pro- 
jected distance from the cluster centre. By assigning stellar 
masses to a single galaxy within each halo, we have also been 
able to test a variety of models for the star formation his- 
tories of cluster galaxies. Novel aspects of this analysis are, 
first, the comparison with stellar absorption linestrengths 
and second, the separate consideration of bulge and disc 
colours, as opposed to global colours. 

In our best-fitting models, the bulge component of 
galaxies is quenched "internally", i.e. its age is not deter- 
mined by infall into the cluster and can be fit with an old 
SSP. The observed colours of the discs of field galaxies re- 
quire that, prior to infall, the star formation rate declines 
with a characteristic e-folding time r prc = 4-5 Gyr. Mod- 
els with no disc quenching fail to explain the redder colours 
of discs near the cluster centre. Models with short quench- 
ing times produce sharp changes in the median colours of 
discs and a stronger cluster-centric radial disc colour gradi- 
ent than is observed. They also overpredict the strength of 
the Balmer lines. Instead, we favour a model in which the 
disc continues to experience exponentially declining star for- 
mation rates with shorter characteristic r post — 2-3 Gyr. 
This is suggestive of a gentler mechanism (e.g. "strangula- 
tion" ) as opposed to rapid ram-pressure stripping of the cold 
gas. 
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